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Abstract 

We study inflationary models that produce a nearly scale-invariant power spectrum while breaking 
scale invariance significantly in the bispectrum. Under most circumstances, such models are 
finely-tuned, as radiative corrections generically induce a larger signal in the power spectrum. 
However, when scale invariance is broken collectively (i.e., it requires more than one coupling to 
break the symmetry), these radiative corrections may be suppressed. We illustrate the features 
and limitations of collective symmetry breaking in the context of resonant non-gaussianity. We 
discuss two examples where oscillatory features can arise predominantly in the bispectrum. 
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1 Introduction 



One of the remarkable features of inflationary cosmology is the robustness of the (nearly) scale 
invariant power spectrum. Despite a vast number of mechanisms that give rise to inflation, their 
predictions for the power spectrum typically differ at the percent level. By contrast, higher Ap- 
point correlation functions can have vastly different behavior and may allow us to distinguish 
these mechanisms. 

The origin of this similarity is an approximate time translation symmetry of the effective 
Lagrangian for the fluctuations. Scale invariance arises in the presence of a time translation 
symmetry because each mode experiences the same history and freezes with nearly the same 
amplitude. This symmetry arises from a shift symmetry for the inflaton, (fi —> eft + c, which 
implies the theory is insensitive to the location of the background field on its potential. 

Shift symmetries in models of inflation have a long history. In the context of slow-roll 
inflation, it is well known that an approximate shift symmetry is required to protect the small 
mass of the inflaton from radiative corrections and from dangerous irrelevant operators (see e.g. 
[1, 2]). Additional work [3-9] has provided explanations for such a symmetry from ultra-violet 
(UV) physics, yielding concrete realizations of these older ideas. 

Recently, it has been suggested that only a discrete shift symmetry is required for viable 
models of inflation [10-13]. In the absence of a continuous shift symmetry, scale invariance may 
be significantly broken by sinusoidal features in power spectrum 1 . In general, one also expects 
large deviations from scale invariance in other N-point correlation functions. Nevertheless, it has 
been shown in the context of single field inflation that the dominant signal of these models is the 
oscillating feature induced in the power spectrum [13]. 

Ongoing observations will constrain (and possibly detect) correlation functions beyond the 
power spectrum, including the bispectrum and trispectrum. When discussing these signatures, 
one typically assumes that all the N-point correlation functions are scale invariant at the same 
percent level accuracy as the power spectrum. This assumption is well motivated by the shift 
symmetry that explains the scale invariance of the power spectrum. If we were to break the shift 
symmetry by the interactions of the inflaton, one would expect radiative corrections to induce 
the same level of breaking in the power spectrum. 

The question we will investigate is whether violations of scale invariance can, in principle, be 
dominated by non-Gaussian correlation functions. We will address this question by constructing 
models where the power spectrum is protected from large radiative corrections by collective 
symmetry breaking 2 . This can be achieved by enlarging the symmetry group of the action such 
that scale invariance is protected by several, independent symmetries. An individual interaction 
may break some of these symmetries, but scale invariance is maintained as long as a subgroup 
is preserved. Therefore, any scale dependent correlation function must be proportional to the 
product of all of the couplings necessary to completely break the symmetry. As a result, these 

1 There is also the interesting possibility that the continuous symmetry is restored in the infrared (IR), despite 
being broken explicitly [14-16]. These examples involve dissipative effects but give rise to scale invariant correlation 
functions. 

2 Collective symmetry breaking is well-known for its use in little-Higgs model building [17, 18]. It has also been 
applied to inflation to explain the absence of large corrections to the inflaton mass [3, 4]. 



2 



effects will occur at the same order in perturbation theory for any N-point function, thus removing 
the preference for the signal to be dominated by the power spectrum. 

In order to realize collective symmetry breaking, we introduce additional fields that trans- 
form linearly (by a phase rotation) in association with the shift of the inflaton. Interactions of 
these fields will explicitly break the continuous shift symmetry, but they will leave an unbroken 
discrete symmetry. As a result, the models we construct will realize a generalization of resonant 
non-Gaussianity [11, 13, 19-21]. These additional fields contribute to the curvature perturbation 
through interactions that convert the isocurvature fluctuation into an adiabatic fluctuation, dur- 
ing inflation. We discuss two classes of models: one where the conversion is perturbative and is 
described by quasi-single field inflation (QSFI) [22] and the other where the mixing is strong and 
is described by a single effective propagating mode. In the two examples we present, there are 
regions of parameter space where the signal is largest in the bispectrum. 

The organization of the paper is as follows: In Section 2, we will review resonant non- 
Gaussianity and its signatures. In Section 3, we show how collective symmetry breaking can 
reduce the resonant signal in the power spectrum, with emphasis on the case of QSFI. In Section 
4, we compute the signals of the perturbative, QSFI example for both the power spectrum and 
bispectrum. In Section 5, we present an example where the mixing with additional fields is strong 
and compute its signatures. Section 6 contains a further discussion of the results. We review the 
relationship between shift symmetries and scale invariance in the Appendix. 

2 Resonant Non-Gaussianity and Radiative Corrections 

In this section, we briefly review resonant non-gaussianity [11] and its origin from a discrete 
shift symmetry in the effective field theory of inflation, following [13]. Resonant non-gaussianity 
occurs when an interaction involving the scalar metric fluctuation, £, is oscillatory in time (e.g. 
Hint ^> A* cos(w^f)C 3 )- When the frequency of this oscillation, w*, is large compared to the Hubble 
scale, H, the correlation function is determined at the time when the modes have energy u ~ to*. 
The time dependence of the oscillation is translated into scale dependence of the correlation 
function through the time of resonance, a(t ies ) ~ k/ca+. The precise relationship between scale 
invariance and time translations is reviewed in the appendix. 

It will be useful to express these results in terms of the Effective Field Theory (EFT) of 
Inflation. The action of the EFT must be diffeomorphism invariant; explicit functions of time 
may appear but will be accompanied by a goldstone boson, tt, that transforms non-linear ly under 
the symmetry [23, 24]. We will work in the decoupling limit (H — > with M^H fixed) such 
that the action is an explicit function of t + tt with ( = —Hit. The goldstone boson transforms 
as tt — > tt — £ (t,x) under a general diffeomorphism x M — > + £^(t, x) which leaves the action 
invariant. 

The Lagrangian for tt is determined at lowest order by the requirement that we have an 
FRW solution to Einstein's equations with tt = 0. It was shown in [23, 24] that the action takes 
the form 

C = -M p 2 ! (3H\t + tt) + ff) + M^Hd^t + vr)^(i + vr) + (t + tt) [0 M (t + ir)&*(t + vr) + l] 2 
+ ^M 3 4 (t + vr) [0 M (t + vr)^(t + tt) + l] 3 + ... . (2.1) 
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If the couplings M^{t + u) are time independent, then the action has a continuous shift symmetry 
under tt — > ir + c. The shift symmetry can be broken to a discrete subgroup if any of the couplings 
are periodic in time. For example, in simplest model of slow-roll inflation, this can be achieved 
using 

Hit) = H ST (t + vr) + H osc (t + tt) sin(w* (t + it)). (2.2) 

In this case, the continuous shift symmetry is broken by H osc to a discrete one tt(x, t) —> tt(x, t) + 
2tt juj*. The discrete shift symmetry allows for self interactions of tt that do not contain derivatives. 
Up to slow-roll corrections, these take the form 



£ n (vr) ~ --^M p 2 ! + 3HH^ vr n . 



Using the above interaction Lagrangian we can calculate the higher correlation functions 



(2.3) 



(2.4) 
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(2.5) 
(2.6) 



We have defined a = cj*/H S> 1 and e osc = °^ 3C , with k+ as the pivot scale. From the quadratic 
action, we find oscillatory imprints on power spectrum, 



7T\l/2 



£osc « 1//2 sin (aln(2fc/fc*)) 



(2.7) 



As it is shown in [13], most of the scale-dependent signal is in the two point function. 

At the classical level, we are free to write the action such that the only oscillatory coupling 
involves more than two powers of £. In this sense, one could imagine generating oscillating non- 
Gaussian correlation functions without altering the power spectrum. However, when radiative 
corrections are included, there is no reason to expect this result to hold. As a concrete example, 
we will consider the interaction H mt = /xcos(u;*i)£ 3 , which arises in the effective theory through 
the operator 



Hi 



hit 



M| cos {u+(t + iv] 
3! 



[^(t + 7r)^(t + vr) + l] J 



M 3 cos(^) [8 . 3 + 12 ^ 7r . 2 + ^ ] 



Mf sin(w*£) 
3! 



8u±tttt 3 + . . .] . 



(2.8) 
(2.9) 



Here, and throughout the paper we will assume a = oj±/H > 1. To compute the radiative 
corrections, we will follow the strategy of [25], where we will estimate the loop contributions 
using a hard cutoff at oj = A. This cutoff is taken to be near the strong coupling scale, A 2 < 
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SL eff . D ^J.Lo ^^ 2 (^(t + 7r))(a^(t + ^)^(t + 7r) + l) 2 ~M2 1 |ij|cos 2 (o;^)^ 2 + ... . 



^■strong = 47r(JW? 1 |fi r |)2/M3, as this is the scale where our EFT is breaking down. Because the 
correlation function is determined at energies of order to*, we further require < A for our 
calculations to be under control. 

We are most interested in the radiative correlations to the quadratic action for tt. The 
largest contribution arises from two insertions of Hi n t that renormalize tt 2 via a 1-loop diagram, 
yielding 

Mj A 4 
(M^\H\) 2 16^" 2 

(2.10) 

We see that we have an order one correction to the kinetic term of tt that depends explicitly on 
time. 

It is straightforward to estimate the ratio of the signal from the resonance in bispectrum 
to the signal in the power spectrum by £3/ £2- This ratio should be evaluated at the resonant 
frequency, oj ~ w*, giving us 

(S/N) 3 A, Mj 

(S/N) 2 ~ £ cfr ~ M 2 | | ?r| w ~ A 2 trong < • (2 - Uj 

Here we have dropped the cos(w*i) since it is of order one (for our purposes). The last inequality 
is essentially the result of [13], which states that the signal is dominated by the power spectrum 
when the effective field theory is weakly coupled. 

3 Collective Symmetry Breaking of Shift Symmetries 

As we saw above, given a generic action with some time dependent coefficient, all operators 
will acquire time dependent couplings of comparable size through radiative corrections. Now we 
would like to understand if there are circumstances where this is not the case. 

Scale invariance is the result of a shift symmetry. This suggests that we can use this 
symmetry to protect the effective action from large radiative corrections. Here we will use the 
idea of collective symmetry breaking to ensure that corrections to the power spectrum arise at 
higher order in the couplings. 

Collective breaking of a shift symmetry arises most naturally in multi-field models of infla- 
tion. We would like tt to shift tt — > tt + c under more than one U(l) such that different U(l)s are 
only distinguished by the transformation(s) of the additional field(s). We can only measure the 
correlation functions of (,/tt, so the transformations of these additional fields are inconsequential 
for scale invariance. To break the shift symmetry, one will need to break all of these U{l)s and 
doing so may require several non-zero couplings. As a result, scale dependence in the tt correlation 
functions must be proportional to all of these couplings, regardless of which TV-point function we 
are computing. 

The most straightforward implementation of this idea arises when tt mixes with an additional 
field, a. For now, let us focus on the case where this mixing is perturbative, namely quasi-single 
field inflation (QSFI) [22]. We still study the observational consequences of these models in 
Section 4 and the case of strong mixing in Section 5. For the remainder of the section, we will 
use QSFI as the primary example only for simplicity. 
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In QSFI, additional light fields contribute to the correlation functions through interactions 
during inflation (rather than through reheating / curvaton effects at late times). To produce a 
signal, we require an additional field a with both self-interactions and a mixing interaction with 
7T. The non-gaussian correlation functions of a give rise to non-gaussian correlations of £ through 
the conversion of a into £ during inflation. The contributions to the various iV-point functions 
are illustrated via the appropriate Feynman diagram in Figure 1. 



a 



\ 

a 1 a 



Figure 1: The leading perturbative contributions to the power spectrum (left) and bispectrum (right) of the 
curvature perturbation £ that arise in quasi-single-field inflation. 

As a concrete example, we will consider the following Lagrangian 

C = M^(3H(t + vr) 2 + H(t + vr)) + M^Hd^t + ^(t + vr) - d^o 

+M pl \H\ l ^ 2 d II (t + -K)d^(pa + pa)[d v {t + K)d v (t + vr) + 1] (3.1) 

Here, a is a complex scalar field, while f3 is a complex, dimensionless coupling and is some fixed 
frequency. For generic values of the couplings, this action explicitly breaks the shift symmetry 
7r — > tv + c to a discrete symmetry tt — > ir + 27rn/w* with n£Z. 

In order to see that this breaking is "collective", first consider the action when /3 = 0. In 
this case, our Lagrangian possesses a exact symmetry, U(l) a where tt — > tt + c and a — > e~ tu> * c ^a. 
Alternatively, we may take (5 ^ but send either M — > oo and / or Wj 0. In this limit, the 
action possesses an exact symmetry, U(l)b under tt — > ir + c with a unchanged. Therefore, the 
shift symmetry of ir (i.e. scale invariance) is only broken in the presence of both w )f 3 ^ and 
finite M 2 . 

The consequences of these symmetries are most apparent if we promote /3 to a background 
field such that /3 — > e - luJ * c / 3 p under U(l) a . As a result, any correlation function of £ containing 
the oscillation will be proportional to ( (3 3 e iu '* tres such that the correlation function is formally in- 
variant under the combined shift and dilatation. Any oscillatory feature in the power spectrum or 
the bispectrum is therefore proportional to /3 3 . This feature ensures that there is no enhancement 
of the signal in the power spectrum for small /5. 

In order to illustrate the implications of the collective breaking, let us consider the renor- 
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C => 3^ L,o,,, lT v M/ j W + h.c] + lTVl1/ ,^ e (/3 2 e^^)d + h.c.)] • (3.3) 



malization of a 2 . Expanding (3.1) in powers of the fields, we find the Lagrangian 

c = -\d^ c d^ c - d^oWa + ^ c {fib + pb) 

where it = \f2M p \\H\ l / 2 TT c and the ... include all terms that will give equal or subdominant 
contributions to the signal to noise. Radiative corrections to the above action will generate the 
following operators 

A \ - [f3e^ t+ ^& 2 + h.c] + — . 

l M 2 {M vl \H\y 2 y J M 2 (M pl \H\V 2 ) 

Both of these operators will contribute an oscillatory signal in the power spectrum. However, 
unlike the previous section, these corrections are suppressed by /3. As explained above, by 
promoting /3 to a background field, we see that operators containing fewer powers of a must 
contain additional powers of (3, in order to maintain the invariance of th.6 action under £/(l) a . 
In writing the second term in 3.3, we have dropped a tadpole for a of the form 

^tadpole ~ -cj^^e^+^a + h.c.) . (3.4) 

Because we are studying the fluctuations, we simply cancel this tadpole by inserting the appro- 
priate counter-term into the bare Lagrangian, <5£ C ounter ~ — ^Aadpole- O ne might worry that the 
presence of the tadpole suggests that a should have a potential with a natural scale set by the 
/i ~ [^uj+jp] 1 / 3 . Due to the suppression by /3 2 , the scale /i can be made to satisfy \i <C H 
without suppressing the signals we will discuss below. 

For the reader who is more familiar with the application of collective symmetry breaking to 
Higgs model building [17, 18], the results in this section may appear slightly unusual. Specifically, 
in a renormalizable theory like the standard model (and extensions thereof), collective symme- 
try breaking is expected to remove power-law divergences. Here, we have suppressed radiative 
corrections by various small parameters, but we have not changed the power of the cutoff. Here 
our interactions are irrelevant and we generally expect power-law renomalization to arise. 

In order to produce a complete model, we still require a UV completion above the strong 
coupling scale that respects the approximate symmetries discussed here. It is also possible that 
such a UV completion would become important at a scale well below the strong coupling scale, in 
effect, lowering the effective cutoff A. Lowering the cutoff would weaken the signal in the power 
spectrum and will be relevant to the generalizations we will discuss in Section 5. For further 
discussion, see [26, 27]. 

4 Perturbative Mixing : Resonance in Quasi-Single Field Inflation 

In the previous section, we showed how collective symmetry breaking can weaken the size of 
radiative corrections to the quadratic action. We are now in a position to compute the actual 
signatures produced by the action 



+ 7 fr[e**- t 6* + h.c] + ^ 2M ^ 1/2M2 [Se^Wd^a + h.c] + £ rad . , (4.1) 
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where £ ra d. are the radiative correlations to the effective action, including those in equation 3.3. 
We first estimate the signal to noise generated in both the bispectrum and power spectrum, 
followed by a detailed calculation of the full momentum dependence and signal to noise. 

4.1 Estimating the Signal to Noise 

The oscillatory features in the power spectrum are controlled by the two operators in (3.3). As 
both will contribute at the same order, let us pick one as an example, namely, 

A 4 

C 2 ^^h --r- W^a 2 + h.c] . (4.2) 

^ M 2 (M pl \H\ 1 / 2 Y J V ' 

Without loss of generality, we may take (5 to be real, such that 

A 4 

/3cos(cj*i)o-! . (4.3) 



32tt 2 



M 2 {M p i\H\y 2 ] 



where ctr = ^(<r + <r). We want to compare this signal with the leading contribution to the 
bispectrum, which arises from 

£ 3 = vW cos M)4 • ( 4 - 4 ) 

These operators contribute to the power spectrum and bispectrum of £ respectively as a result 
of the mixing term f3n c &R, which convert the ctr fluctuations into fluctuations of £ = —Hit. 

We will estimate the ratio of the signal to noise in the bispectrum, (S/N)s, to signal-to- 
noise in the power spectrum (S/N)^ '. In a typical single field model, we would estimate this by 
(£3) / '(£<r c ') evaluated at horizon crossing, go ~ H . There are two important subtleties that arise 
in these models. First, in the case of resonance, the signal is generated at the resonant frequency 
00 = uj+ rather than horizon crossing. Second, we want to estimate the signal in Q correlation 
functions, not a correlation functions. The mixing interaction is responsible for converting a a 
to a C with strength 3 /3a, where a = 00+/ H. Therefore, our ratio of the signals-to- noise is given 
by 

(S/N) 3 a 3 /? 3 £ 3 U^ a{M pl \H\ l l 2 )a R (M pl \H j 1 / 2 )^ 3 
(S/N)f c - ~ a 2 /3 2 £!f c -U = ^ ~ A 4 " K A H ' 1 ' ' 

In order for this estimate to be reliable we require the resonance to occur below the strong 
coupling scale A* ~ M, or oo± < A*. Even if we take the UV cutoff to be the strong coupling 
scale A ~ A* ~ M, the non-gaussian signature dominates when 

M[ — ffg— ) <«*<M. (4.6) 



(M 2 |ij|)V2 

Therefore, if M 2 H 2 < M 2 X \H\ (or M 2 < eM 2 j) there is a range of oo± where the resonant signal 
is mostly non-gaussian. In general A < A*, which would increase the range of parameters where 
the bispectrum gives the largest signal. 



3 Naively, the interaction strength is given by /3. The enhancement by a represents the extra time over which 
the mixing operates relative to conventional freeze-out. This will be shown in detail in the next section. 
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4.2 Calculating the Bispectrum 

In this section we will present calculations of the bispectrum and power spectrum. The calcula- 
tions are nearly identical, therefore we will only show the bispectrum calculation in detail. 

Cosmological correlation functions are given as equal time correlation function in the in-in 
formalism. Following the standard presentation (see [28] for details), we compute correlations 
function in the interaction picture as 



{Texp[i[° flirt (r)dr] (Cint(^i, r ). . Cint(4, r ))rexp[-f f ° fl in t(r)dr]) , (4.7) 

J-oo(l+ie) V ' J—oa{l—ie) 



where we are working in conformal time r. In practice, it is easier to perform this calculation 
after Wick rotation. From the ie prescription, it should be clear that we must rotate the time 
and anti-time ordered exponentials in the opposite way such that we are calculating 

(Texpli/ 1 H- mt (T)dT] (Cint(£i, r ). . Cint(fen, 7b))rexp[-i / L flint (r)dr]) . (4.8) 

However, this is nothing other than the anti-time ordered correlation function in Euclidean time 

/I OO "J" TO 
flint (r)dr]\) (4.9) 
-joo+ro 

_ / -, -, roo—vnt . 

= (TUmt(ki,To)..Cmt(kn,ro)exp[- H mt {-iT E )dT E ]\) . (4.10) 

V J -oo+ro 

where r = ite- It is now straightforward to calculate the bispectrum 4 , provided we use the 
anti-time ordered, Euclidean Green's function 

{f(a{T E ,k)(T{r' E M))) = ^(l + h\T E - r' E \ -klT E T'E)e- k ^- T, E\(2TT) 2 5{k l + k 2 ) . (4.11) 

The Green's functions for £ are the same, up to the overall normalization which we define such 
that C(r -> 0) = ^ with A c = 2.2 x 10" 4 [29]. 

We will compute the bispectrum, setting tq = for simplicity and writing the Euclidean 
conformal time te — > r. Our interaction Hamiltonian is given by 

fl int = (flr^Mpilflf/ 2 ™ + cos(Lu*t)d 3 } , (4.12) 

where t = — fl _1 [log(Hr) — |] in de Sitter (A = 0) and where have dropped the index on ctr — > a. 
The bispectrum is given by 



3! COS(W 4 ) . . . , x lx 

Tl)] [ 75 (flr 4 ) 4 M 2CT -^' 7 -^ CJ -fc3 ^ 

(4.13) 



4 We thank Kendrick Smith for this suggestion. 
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Here we have implicitly contracted the fields, integrated over the momentum and dropped the 
overall momentum conserving delta function (27r) 3 <5(/ci + k 2 + k^). Using the anti-time ordered 
Green's functions, this becomes 

A 3 r°° B~2_2 



2 ° V3/2, 3/2, 3/2 

nji fto Ajq 



/CO f[ T f°° 
dn— ± cos(alog(#T 4 ) - ian/2) \[ / dn((3k ie - k ^~ k ^-^) 
-OO • -i J — oo 



1 ^2 ^3 ""^ i=l 



= |( fe fc C fc ^ ) / dr 4 r| cos(alog(^r 4 ) - joTT^e"*™ + ^|r 4 |) . (4.14) 

This integral is straightforward to evaluate in general, although we are primarily interested in the 
leading behavior in a > 1. Because the dominant contribution to the integral is from r ~ ce/K, 
the leading contribution is from 

1 - K A C1^2-) J ^4|r 4 | 5 e-^M cos(alog(Fr 4 ) - iavr/2) (4.15) 

= |||^ 2 [|(r(6 + i«)^ + r(6- i a)K + -)ef" (4.16) 
+|(r(6 - ia)K +ia + T(6 + ia)K- ia )e~l a \ . (4.17) 
Using Sterling's approximation, up to a ^-independent phase, we find 

(C^CfeC^) ^ |(||^)v / 2^a 11/2 cos(alogK/^)(2vr) 3 5(fc 1 + k 2 + fc 3 ) . (4.18) 

The power spectrum calculation is essentially identical, except that now 

= (Ht)- a [PM p1 \H\ 1 / 2 7t& + 7 cos(c^)a 2 ] , (4.19) 

where 7 ~ 'm^m^\h\^) ^ when it is generated from radiative corrections. Following the steps 
in the previous section, we find for a S> 1 

A 2 

(C^k 2 ) - (^ 2 l)V2^(2)- 3 a^\os(alogk 1 /k,)(2nf5(k 1 + k 2 ) . (4.20) 

Comparing equations 4.18 and 4.20, we find good agreement with our parametric estimate in 
equation 4.5. 

It is straightforward to calculate the signal to noise, following the appendix of [13], including 
all numerical factors and the full momentum dependence. To get an accurate estimate, we 
compute the bispectrum at all orders in a using equation 4.14. This ratio is shown in Figure 2 
in terms of -r and 

The signal to noise for much of the parameter space is surprisingly small given the parametric 
estimates in the previous subsection. This is due to a small numerical coefficient in (-^ ((C 3 )))- 
The suppression is most easily understood by considering the bispectrum 4.18 in the equilateral 
configuration k\ = k 2 = k%. Because K = 3fci in this configuration, the amplitude is suppressed 
by 3~ 6 ~ 10" 3 . For this reason, the signal to noise ratio is on the order of 10 — 100, rather than 
the 10 4-5 that we expected. Nevertheless, the bispectrum dominates over an appreciable range 
of parameters where the theory is under control. 
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Figure 2: Contour Plot for )j f )', /v. . To ensure reliable results for the small values of a, the signal to noise 

\w\ °\* ))) 

includes the full a dependence of the bispectrum given by equation 4.14 (and similarly for the power spectrum). 
We have enforced a > 2 to ensure that the our calculations and the theory are under control. 



5 Strong Mixing : Resonance with a Linear Kinetic Term 

In the previous sections, we looked at one particular realization of the collective breaking. In this 
section, we will discuss a related example to show that the results were not particularly sensitive 
to being in the QSFI regime. 

We will change the model from the previous section by altering the form of the mixing term. 
The Lagrangian we will consider is instead 

£ = M^(3H(t + tt) 2 + H(t + tt)) + M^Hd^t + Tr)d^{t + vr) - d^a (5.1) 
+M pl \H\ l / 2 (pa + pa)[d^t + Tx)d^{t + vr) + 1] 
+ M2 [e^^id^t + vr)c"V) 3 + h.c] . 

Here p is a complex coupling with dimensions of energy. The only change to the above action 
from equation 3.1 is in the second line. Keeping only the necessary operators and canonically 
normalizing tt ^ 2M 2 ^\H\ = tt c , we arrive at 

£ = \d^ c d^ c + d^a + ^c(pa + pa) + ^ [e^a 3 + eT^rj 3 ] + . . . . (5.2) 

When p <^ H, we expect the results to be identical to the previous section with (3 <C 1 after 
replacing /3 — > p/H. Here, we will be interested in the regime where /)> H. 

The dynamics of this model depend significantly on the energy scale [26]. At high energies, 
co 3> p, this is a two field model with tt and a acting as the independent degrees of freedom. 
When uj <C p, the mixing term dominates: p-fra ~ pujiro S> 7i 2 ~ uj 2 tt 2 . At these energies, a is 
not an independent degree of freedom, but rather pa ~ P n is the canonical momentum of ir c . We 
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will be interested in uj* < p such that we are effectively dealing with a single field at resonance. 
In this case, pitcGR is effectively the kinetic term for both ir and a and we may drop tt 2 and a 2 
as they become irrelevant operators. Therefore, the effective action is given by 

C = p-k c a R - i<9j7r c dV c - \di<j R d l a R + ^J^ 2 cos(ujj)a R + ... . (5.3) 

where we have taken p to be real and a R is the real part of a. The imaginary part of a will not 
play role so we will again drop the subscript, a R — > a. 
The quadratic equations of motion are given by 

k 2 

pn k(t = , (5.4) 

cr 

—pa + 3Ha r ^tt = . (5-5) 

a 1 

The first equation is just a constraint that requires a = a 2 p-k c /k 2 , so we see again that a is not 
an independent degree of freedom. Plugging this back into the equations, one can solve for the 
the mode function of tt c and are given by [26] 

^)=^ 2 ^(-) i/2 ^;i(if(^) 2 ) ) (5.6) 

such that 

ir c (k,T) = a^ + alut . (5.7) 
In terms of the parameters of the model, the power spectrum (for M — > oo) is given by 

where ' indicates we have removed the (27r) 3 <5(/c + k'). 

The radiative corrections discussed in Section 3 are identical provided we substitute j3 — > 
p/A. However, the model is not free from dangerous correction, even when the shift symmetry 
for 7r is exact. This is because the action no longer possesses a shift symmetry for a and we must 
worry that a large mass m 2 a 2 will be generated. We require m <C p for the above dynamics 
to hold. Fortunately, due to the modified dynamics at low energies, this mass can be made 
sufficiently small without modifying the model below the scale p [27]. 

5.1 Estimating the Signal to Noise 

Before doing a detailed calculation, we will estimate the signal in the bispectrum compared 
to the power spectrum. The power spectrum is scale invariant at leading order, but receives 
contributions from radiative corrections. The largest effect comes from the renormalization of 
the operator 

Cf c - = 3 ^ cos^d 2 . (5.9) 
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The largest contribution to the bispectrum will be generated by £3 = cos(cj^t)<7 3 . As 

before, we can estimate the ratio of the signal to noise to the bispectrum to the power spectrum 



as 

1 A3 



(S/N) 3 £3 , . 7^P & , MpiH 1 / 2 a 



(S/N)f c - _ ^_ & 2 lu - u * A 3 p l 

M vl H l l 2 u 2 tt c 



A 3 (k 2 /a 2 ) 



(5.10) 



(5.11) 



To complete the estimate, we need to understand the behavior of ir c . Because of the non- 
relativistic kinetic term, this is not really the "canonically normalized" field in the sense that its 
mode functions are not purely functions of energy and momentum. Instead, tt c = TTcp 1 ^ 2 is the 
canonical field, but is has units [tt c ] = [^] 3//2 . Using this observation, with u = k 2 /(a 2 p), we find 

(S/N) 3 MpigVg co 2 jVyjV 2 p /^n7/4 

(S/N)f c - ~ A 2 p i/2(^)i/2 |w -^ ~ A 2 AV/9/ ' 1 j 

The bispectrum will produce the largest signal provided that 

A 2 A 1 



Mpitf 1 / 2 p 



7 < ^ < p (5.13) 



The upper bound w + < p < A is an assumption of the setup. Therefore, a necessary (but not 
sufficient) condition for satisfying (5.13) is that A 2 < M V \\H\ 1 / 2 . 

5.2 Calculating the Bispectrum 

In order to determine the bispectrum to leading order, we must calculate 

where d 3 (r) = / ^ s ^^^^&^ 1 a^ 2 a^(2iT) 3 5 3 (pi + p 2 +^3) and a = uj+/H. Because a = 
pa 2 7r c /k 2 , we can determine & as an operator by differentiating 7r c , using (5.7). Therefore, we 
need to evaluate 

where a % = a % a % + a\at. 

The integral in (5.14) is dominated by the saddle point at r 2 = j^jj with K 2 = kf + k^ + k 2 . 
Because we are interested in the behavior when a ^> 1, we can keep only the term with the most 
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factors of r, as it gives the largest power of a. The leading contribution to (5.14) is given by 

3! H^k^h) 3 / 2 tt^A 3 



y/2M 2 p 6 
x2Re / drr 19/2 cos(alog|^r|)]J^ 

J— oo(l-ie) j =1 

3!# 2 H^ik^h) 1 / 2 A3 



( 2 ) fl^fi-V) 

3/4^2 p ^* T J 



9 9/2 



Jl/2 



V2M 2 p 

/•o _4# 
x2Re / dTT 13/2 cos(a log | Ht | )e * 2 P 

J — oo(l— ie) 



(C^C^C^) - p6 83/2 (5-16) 

/■o JL , , 

(5-17) 

(5.18) 
(5.19) 
(5.20) 

., - _\ y /I ( ! V UfM f. ! LV l\ i h + ! . (5.21) 

M 2 4^3/4^15/2 C V / *7 ' V 7 

where we have absorbed a fc-independent phase into the definition of the pivot scale fc*. 

The oscillatory signal in the power spectrum follows from a nearly identical calculation. 
Following the same steps as above, we find 

- ( ~ , ^1 ; - )2Ref ^el^fcf r[2 - if]) (5.22) 



(Kr) 2 



15 .Q- 



M 2 2 1 /4 /0 3/4 i ^15/2 C V L 

-A^V 71 " cos(alog K/K+) , 



32 " A: 3 VMp^l^M 2 / V4p 
A 2 , HA 3 _ , 3/2 



32 - 2 kf \M pl \H\V 2 M 2 



7r(a) ' cos(a log ki/k+ 



(5.23) 



1 1 < 



o.o - 




0.00 



Figure 3: Contour Plot for 



(f «c 3 ))) 



with q = uj* I H = 10. The signal to noise from the bispectrum and the 



(#(*«=») 

power spectrum was calculated using the full a dependence from equations 5.14 and 5.15. 

As in the previous sections, we can use the results of [13] to give a more precise estimate of 
the signal-to- noise. To achieve accurate results for moderate values of a, we determine the full 
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a dependence by including both terms from (5.15) when computing the bispectrum and power 
spectrum. The signal to noise ratio is shown in Figure 3 as a function of ^ and ? for fixed a = 10. 
From our naive estimates, we might have expected a smaller signal than what we found the QSFI 
example in the previous section. Specifically, the signal in the bispectrum is only enhanced by 
a 7 / 4 relative to the power spectrum compared to a 3 in QSFI; yet, the signal to noise estimate gives 
comparable results for physical values of a. The main reason is that numerical suppression of the 
signal to noise is smaller than in the previous section. We can again estimate the suppression by 
looking at the equilateral limit k\ = k% = h^. In this configuration, K = \Jk\ + fc| + ~ \/3fei 
and our bispectrum is suppressed by 3 -15 / 4 ~ 10~ 2 . This is larger by a factor of ten from the 
QSFI example and explains why the results are comparable for a ~ O(10). 

6 Discussion of results 

We studied the possibility that collective breaking of a shift symmetry could give rise to a large 
oscillatory signal in the bispectrum without introducing larger oscillations in the power spectrum. 
We constructed two examples where the radiatively induced signal was suppressed relative to the 
bispectrum. The common feature of these models is the presence of an additional scalar field 
that transforms by a phase in association with the shift of the inflaton. 

Scale invariance in the power spectrum and the bispectrum is broken in our models through 
oscillations in momentum space that arise as a result of resonance [11, 13]. Due to the rapid 
oscillations, it is unclear that the signal to noise in the CMB is well approximated by the signal 
in the primordial correlation functions. Given that the signal in the power spectrum is smaller 
than in the bispectrum by a factor of ten to a hundred, a more careful analysis of the signal in 
the CMB is required. 

The broader motivation for this work was the question of whether non-gaussian correction 
functions can naturally contain generic, order one violations of scale invariance (i.e. including 
non-oscillatory shapes). This is relevant to future observations of scale dependent bias [30, 31] 
and ^-distortion [32], where there is a degeneracy between the effects of massive fields [22, 27] or 
excited states [33, 34] and a scale dependent amplitude for /^£ al [35, 36]. It was crucial to our 
above constructions that we only broke the shift symmetry to a discrete subgroup. Power law 
violation of scale invariance would require that the shift symmetry is completely broken. The 
approach of collective symmetry breaking used in this paper would not protect such models from 
large radiative corrections. It remains an open question whether large, power-law violations of 
scale invariance can arise in the bispectrum without fine-tuning. 
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A Scale Invariance and Shift Symmetries 



In this appendix, we will briefly review the relation between shift symmetries and scale invariance. 
An FRW universe possesses a dilatation symmetry that arises from a large diffeomorphism. This 
symmetry is realized non-linearly by the curvature perturbation, £, and takes the form in co- 
moving coordinates [37] (see [38-40] for further discussion) 

x->(l + A)x C -> C + A(l + x -d s () , (A.l) 

where A is an infinitesimal parameter. Due to the presence of this symmetry, there must exist 
a solution (,{x,t) = Co> which is the limit of a finite momentum solution [37]. This ensures 
that correlation functions of long wavelength modes are time independent. However, because 
the transformation acts non-linearly on £, the correlation functions are not required to be scale 
invariant. 

Formally, the action for £ is also invariant under the rescaling x — > Xx and a — > X~ 1 a. 
However, the scale factor, a(t), is itself a function of time and therefore this symmetry is broken 
by any time dependent coupling. For example, in de Sitter, by writing t = H~ 1 loga, we see 
that a time dependent coupling X(t) — > A(a) breaks this symmetry. On the other hand, if the 
action is invariant under t — ¥ t + c, then our rescaling symmetry is also preserved. Under the 
rescaling, Cfc(t) — > A~ 3 Cfc(i). In addition, must be time independent as k — > which implies 
Cfc(^) ~~ ^ 3 m this limit. In other words, the correlation functions of £ are scale invariant in the 
long wavelength limit, if the action possesses a time translation symmetry. 

The physics of freeze-out also makes this connection clear. The correlation functions are 
typically determined by the couplings in the effective action for £ at the time of horizon crossing, 
e.g. for slow roll when k = a(i*)-ff(i*). When we observe modes at different values of k, they 
were sensitive to couplings at different times t*(k). Hence, the time evolution of the couplings is 
translated into scale dependence through the freeze-out time t*(k). In resonant non-gaussianity, 
the role of freeze-out is played by the time when the mode crosses a resonance, k = u>*a(i res .). The 
general conclusion will be the same in both cases, couplings that are time dependent during in- 
flation introduce scale dependence in late time observables through the freezing of the correlation 
functions at some fixed physical scale. 

Although it may not seem clear from this presentation, this is closely related to having 
shift symmetry for the inflaton. This can be most easily understood in the context of the effec- 
tive theory of inflation. During inflation, the time diffeomorphisms are broken by the inflaton 
background, and are only realized non-linearly in the theory of the fluctuations. The time dif- 
feomorphisms are maintained through a shift of the goldstone boson, it, provided that the action 
is a generic function of t + ir. As a result, if the action is time independent, then the symmetry 
t — > t + c implies there is a symmetry under which (t + n) —> (t + 7r) + c. This symmetry can also 
be interpreted as the statement that the action is invariant under a shift symmetry tt — > it + c. 
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